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ABSTRACT 



Of 



The formation mechanism of neutron stars with extremely large magnetic field 
strengths (magnetars) remains unclear. Some formation scenarios predict that 
magnetars should be born with extremely high space velocities, > 1000 km s -1 . 
Using the Long Baseline Array in Australia, we have measured the proper motion 
of the intermittently radio-bright magnetar PSR J1550-5418 (IE 1547.0-5408): 
H = 9.2 ± 0.6masyr _1 . For a likely distance of 6 ± 2kpc, the implied transverse 
velocity is 280^} 3 ; okms~ 1 after correcting for Galactic rotation. Along with the 
f=s 200 km s -1 transverse velocity measured for the magnetar XTE J1810-197, 
this result suggests that formation pathways producing large magnetic fields do 
not require very large birth kicks. 
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Introduction 



Magnetars, originally identified via their bright and variable high-e nergy emission, are 
neutr on stars largely powered by the decay of ultra-strong magnetic fields (IDuncan fc Thompson 
1992 ). Known examples of magnetars^ - both anomalous X-ray pulsars (AXPs) and soft 
gamma-ray repeaters (SGRs) - are young, but rotate s lowly because their st r ong d ipole 
fields (B x ~ 1Q 14 ~ 15 G) hay e spun them down quickly. IDuncan fc Thompson! (119921 ) and 



Thompson fc Duncan! (11993! ) suggested that such large fields can only be created through 
dynamo action in proto-neutron stars with ~ 1ms spin periods, which is shorter than the 
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convective overturn time. Searches for evidence of such fast initial rotation via its enhanced 
contribution, up to 10 52 erg, to the explosion ener gy of the supernova remnants (SNRs) host- 
ing m agnet ars, have generally come up negative fjVink fc Kuiperll2006l ; lHalpern &: Gottheli 
20101 ) . On the other hand, it is possible that the diffusion of such rotational energy through 
the envelope of a supernova can exp lain some of the most lu minous SNe light curves that 
have been observed in recent years (IKasen fc Bildsteinl 120101 ) . The birth of a millisecond 
magnetar is also considered a possible channel for gamma-ray bursts. 

A possible consequence of the large internal magnetic fields (B > 10 15 G) or fast initial 
rotation (P < 1 ms) of magnetars is a very high kick velocity. Processes that are ineffective 
in ordinary neutron stars, but may produce v > 1000 km s -1 recoils in magnetars, include 
anisotropic mass l oss in a wind or jet, the electromagnetic rocket effect, and anisotropic 
neutrino emission ( IDuncan &: Thompson! Il992l ; iThompson fc Duncan! Il993l ) . On the other 
hand, the central location of some m agnetars within SNR s has led to the inference that as 
a class their velocity is < 500 kms -1 ( jGaensler et al.ll200ll ). closer to that of the population 
of ordin ary pulsars. The superposi tion of SGRs on massive star clusters or giant molecular 
clouds ( jCorbel fc Eikenberryl I2004J ) that are considered their likely birth places also offers 
circumstantial evidence that their velocities are unremarkable. 

Attempts to actually measure velocities via X-ray or infrared imaging remain unful- 
filled, and only the first radio-detect ed magnetar, XTE J1 810-197, has a reported proper 
motion. Using VLBA observations, iHelfand et al.l ( 120071 ) obtained a transverse velocity 
of m 200 kms" 1 for XTE J1810-197, slightly below the average for young neutron stars. 
Although the unknown radial component of velocity makes it impossible to draw firm con- 
clusions, this result is not supportive of theories of magnetar formation mechanisms that 
require a high initial velocity. However, more measurements are required in order to actually 
measure the initial velocity distribution of this exotic class of neutron stars. 



The discovery of a second radio-emitting magnetar, PSR J1550-5418 (jCamilo et al. 



20071 ). provided the opportunity for a second proper motion measurement. Here we present 
Long Baseline Array observations and the resulting proper motion for this 2 s pulsar (also 
known as IE 1547.0-5408). Our VLBI observations and results are presented in Sections [2] 
and El and the implications for the birth velocity of PSR J1550-5418 and magnetars in 
general are discussed in Section HJ 
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2. VLBI observations 



We observed PSR J1550-5418 five times over a period of 2.5 years using Australia's Long 
Baseline Array (LB A). The LBA consists of six antennas: the Australia Telescope Compact 
Array (ATCA), Parkes, Mopra and Tidbinbilla antennas operated by CSIRO Astronomy and 
Space Sciences, and the Hobart and Ceduna antennas operated by the University of Tasma- 
nia. The first observation, obtained using a target of opportunity trigger in 2007 August, was 
made at 2.3 GHz, and revealed that PSR J1550-5418 suffers from severe scatter broadening 
due to the interven ing; ionized interstellar medium (already measured in the time domain by 



Camilo et al.l 120071 ) . PSR J1550-5418 was detected on the shortest baselines only, and the 



scatter broadening was estimated to be approximately 80 milliarcseco nds (mas) at 2.3 GHz . 



This scatter-broadening exceeds the prediction of the NE2001 model (ICordes fc Lazidl2002l ) 
for the location of PSR J1550-5418 by a factor of 5. Accordingly, the remaining observations 
were made at 8.4 GHz, where the effects of scatter-broadening are much less severe and the 
astrometric precision much higher. 

Observing setups varied between epochs due to resource constraints and equipment 
upgrades, with observation durations between 5 and 12 hours, and total bandwidths of 
32 MHz or 64 MHz (dual polarization). Table [1] summarizes the 8.4 GHz observations, which 
were used for the astrometric results detailed below. 

At the commencement of observations, the sparse calibrator grid in the southern skies 
meant that no known phase reference sources were available within 6° of PSR J1550-5418. 
Instead, the observations were phase referenced to the bright and compact radio source 
J1515-5559 (~ 1 Jy at 8 GHz). J1515-5559 was identified as a probable calibrator based 
on lower-resolution data and confirmed as such in the first astrometric epoch. However, its 
relatively large angular separation from PSR J1550-5418, 5.°3, means that the quality of 
phase referencing is much poorer than with other equivalent LBA observations at 8.4 GHz 



Table 1. Summary of astrometric observations 



Observation date 


Participating telcscopcs a 


Observing bands b 


Duration 


Detection S/N 


(yymmdd/MJD) 




(MHz, dual polarization) 


(hours) 




071111/54415 


AT, CD, HO, MP, PA 


8409-8425,8425-8441 


5 


6 


080327/54552 


AT, CD, HO, MP, PA 


8409-8425,8425-8441, [8441-8457,8457-8473] 


12 


11 


090226/54888 


AT, CD, HO, MP, PA 


8409-8425,8425-8441, [8441-8457,8457-8473] 


10 


20 


091211/55176 


AT, CD, HO, MP, PA, TI 


8409-8425,8425-8441, [8441-8457,8457-8473] 


8 


9 



a AT = Australia Telescope Compact Array, CD = Ceduna, HO = Hobart, MP = Mopra, PA = Parkes, TI = Tidbinbilla. 
b Additional observing bands listed in square brackets were only present at AT, MP, PA. 
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(e.g., iDeller et al.l 120081 ). A 5 minute target /calibrator cycle was employed, with 2 minute 
calibrator scans and 3 minute target scans. 

After the first epoch, scans to check phase referencing were included on other known 
calibrator sources within 6° of J1515-5559, and the results used to refine its position. As the 
position for the calibrator source J1515-5559 was refined in successive observations, different 
positions were used in correlation at each epoch, but in the final analysis all epochs were 
corrected to the more accurate position ultimately available in the rfc_2011d VLBI source 
catalog^, (J2000.0) R.A.=15 h 15 m 12. s 67312, decl.=-55°59'32. // 8361, with errors of 5.4mas in 
R.A. and 2.0 mas in decl. This position was derived through observations in the LBA Calibra- 
tor S urvey, a dedicated southern hemisphere campaign for absolute astrometry (IPetrov et al. 
201lh . 



The data were correlated with the DiFX software correlator (IDeller et al.l 120071 ) . using 
matched filtering (gating) on pulse profiles. Pulsar ephemerides were supplied at each epoch 
from accompanying timing observations done with the Parkes telescope, and were used to 
set gates of widths that ranged between 7.5% and 20% of the pulse period. The applied gate 
width depended on the (varying) pulse profile at each epoch and the uncertainty in the pulse 
phase prediction caused by unstable rotation. Due to the relatively high observing frequency 
of 8.4 GHz, no ionospheric correction was applied after correlation. During some of the ob- 
servations, the pulsar flux density varied significantly on tim escales of hours (these variations 
are intrinsic to the magnetar; see lCamilo et al.l 120071 120081 ) . The a mplitudes were flatte ned 
(and the d ata reweight e d accor dingly) using a ParselTongue script (IKettenis et al.ll2006l ) de- 
scribed in IDeller et al.l ( 12009bl ). Correction of the amplitudes in this manner minimizes the 
flux scattered away from the true pulsar position and maximizes the significance of the final 
detection, at the expense of an elongated beam shape due to the time-variable weighting. 
After imaging using natural weighting, the AIPS task JMFIT was used to estimate positions 
and errors. The brightness of PSR J1550-5418 varied considerably between epochs, but it 
was detected in all epochs with a significance ranging from 6 a to 20 a (Tabled]), leading to 
positional fits with a nominal accuracy of ~ 0.2 — 0.6 mas. 

At 8.4 GHz, the deconvolved, scatter-broadened angular diameter of PSR J1550-5418 
was found to be ~ 6.5 mas. In the four astrometric epochs, fits to the semi-major axis 
ranged between 5.5 and 8 mas. This angular broadening leaves very little flux on baselines 
with length > 1000 km, meaning these baselines contribute very little to the position fit. 
Consequently, the longest LBA baselines (to Hobart and Ceduna) contributed relatively 
little to the astrometric accuracy in this project. 



2 http://astrogeo.org/vlbi/solutions/rfc_2011d/ 
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3. Astrometric results 

Due to the large angular separation between the target and the calibrator, the astro- 
metric errors are dominated not by the signal-to-noise ratio but by systematic contribu- 
tions. The predominant systematic effects are introduced by the calibrator position uncer- 
tainty (~ 5 mas), the station position uncertainties (several centimeters) and the unmodeled 
atmosphere (delay equivalent to a path length of several centimeters). Both simulations 
( iPradel et al.l 120061 ) and extrapolation from past results with the LBA ( iDeller et al.l 120081 ) 
suggest that under these conditions, the systematic contribution to astrometric error should 
be of order 1 mas. 

For the astro metric fit, the mag netar parallax was fixed at 0.11 mas, corresponding to a 
distance of 9 kpc (jCamilo et al.ll2007l ). Varying the distance over a wide range of values, from 
2 kpc to infinity, made an insignificant (~ 0.1 a) change to the proper motion fit. Without the 
addition of any systematic error estimate, the reduced x 2 OI t ne fit was 8.0, indicating that 
as expected the formal errors in fitted position considerably overestimate the astrometric 
accuracy. The addition of 1.2 mas (R.A.) and 0.6 mas (decl.) errors in quadrature was suffi- 
cient to obtain a reduced x 2 OI 1-0- The optimal values for the systematic error (minimizing 
the squared sum of the R.A. an d decl. components) were found via iterative minimization, 
following the technique used by iDeller et al.l (J2009aj) . The result of the final fit is shown in 
Table |2J The absolute positional accuracy of the reference position for PSR J1550-5418 is 
dominated by the uncertainty in the absolute calibrator position. The observing history of 
the calibrator source is not long enough to determine if it exhibits time-variable apparent 
position changes (which would corrupt the target proper motion), but in this case the in- 
fluence of such effec ts can be neglected, sin ce the magnitude of such changes are typically 



<0.1masyr" 1 (e.g.. MacMillan fc Mall2007h . 



Figured] shows the positional measurements and fitted motion for PSR J1550-5418. The 
accuracy is considerably better in decl. than in R.A. due to the predominantly north-south 

Table 2. VLBI results for PSR J1550-5418 



Parameter 


Fitted value and error 




R.A. (J2000) a 


15 h 50 m 54? 12386 ± 0. H 00005 ± Of 


00064 


Decl. (J2000) a 


-54° 18' 24." 1141 ± 0."0003 ± 0.' 


'0020 


Epoch of position (MJD) 


54795.0 




Proper motion in R.A., ft a cos S 


4.8 ± 0.5 masyr" 1 




Proper motion in Decl., fA$ 


-7.9 ± 0.3masyr _1 




Total proper motion, ti 


9.2 ± 0.6masyr _1 





Note. — All errors listed arc 1 a confidence level intervals. 

a Cclestial coordinates are given i fit errors ± calibrator position errors. 
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arrangement of antennas in the LBA. 



4. Velocity of PSR J1550-5418 and of magnetars 



Our results provide a secure measurement of the proper motion of PSR J1550-5418, but 
converting this to a transverse velocity requires an estimate of the magnetar's distance. The 
free electron column density obtain ed from the mea s ured d ispersion of the radio pulses, along 
with the electron den sity model of Cordes &: Laziol ( 120021 ). suggests that PSR J1550-5418 is 
located at d = 9kpc (jCamilo et al.ll2007l ). The uncertainty on this estimate is unknown, but 
could approach 50%. An entirely different method relies on modeling three dust-scattered 
X-ray rings observed following a huge outburst of the magnetar that took place in 2009 
Januar y (one month before the third LBA epoch listed in Table [1]). Using the scattered 
X-rays, iTiengo et al.l (120101 ) obtain a best-fit distance of 4kpc, with anything in the range 
4-8 kpc plausible. 

We therefore consider that a reasonable distance estimate for PSR J1550-5418 is d = 
6 ± 2kpc. Along with our measured proper motion of fj, — 9.2 ± 0.6masyr _1 (Table |2J), the 
implied transverse velocity is v± = 260 ± 90kms~ 1 . Correction for peculiar solar motion 
and Galactic rotation using a flat rotation curve and the current IAU recommended rotation 
constants (R = 8.5 kpc, 6 = 220 kms -1 ) adjust s this slightly to 2 80^90 km s -1 , typical 
for the general population of ordinary pulsars (e.g.. lHobbs et al.ll2005l ). This velocity is also 
comparable to th at obtained for the magnetar XTE J1810-197 (v± = 212 ± 35kms~ 1 for 
d = 3.5 ± 0.5 kpc; iHelfand et al.l 120071 ) . The small angular velocity of PSR J1550-5418 also 
implies that an association between the pulsar and th e possible SNR shell G327.24-0.13 



that apparently surrounds it (IGelfand fc Gaenslerll2007l ) rem ains plausible (see ICamilo et al. 
2007 ). Since the characteristic age of the pulsar is < 1400 yr ( Camilo et al. 20081 ) . the pulsar 
could have moved at most ~13" since its birth. This is consistent with its location near the 
centre of the apparent SNR, which has a diameter of 4'. 

With these two proper motions in hand, it is already very unlikely that magnetars as a 
class have exceptionally high space velocities. For example, if we assume that all magnetars 
have v = 1000 kms -1 , then the probability that one will be observed with v± < 212 kms -1 
is 0.023, while the probability that v± < 280 kms -1 is 0.04. Jointly, the probability of both 
observations is 9 x 10 -4 under this assumption. If all magnetars have v = 500 kms -1 , then 
the corresponding individual probabilities become 0.094 and 0.171, with 0.016 for the joint 
probability. We therefore consider v < 500 kms -1 as a reasonable (~ 2.2c) upper limit on 
the typical magnetar velocity. Additional measurements or constraints on magnetar proper 
motions would clearly be desirable to further constrain their distribution of space velocities. 
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Unfortunately, X-ra y observations with Chandra have yielded only upper limits that 



are not very restrictive. Kaplan et al.l (120091 ) find for IE 2259+586 fi < 65masyr l ,v± < 



930 km s 1 assuming d = 3kpc, which is not as small a limit as has already been inferred 



from t he pulsar's central location in its host SNR CTB 109. For SGR 1900+1 4. Kaplan et al. 



(2009) find /i < 54masyr 1 ,v± < 1300 km s 1 assuming d = 5kpc (see also iDe Luca et al. 
20091 ). This result is not constrain ing if an associa t ion of SGR 1900+ 14 with a young star 
cluster at d= 12.5 kpc is accepted JVrba et alihoool : bavies et alJboogh . since SGR 1900+14 



remains within 15" of the center of the cluster . Similarly, the ap parent association of SGR 
1806-20 with a cluster of giant massive stars (IFuchs et al.lll999l ) requires that magnetar to 
have a small velocity. 



The only other known radio-emitting magnetar, PSR J1622-4950 (ILevin et al.ll2010l ). 
is an attractive target for the LBA. But in the meantime, our LBA measurement of the 
proper motion of PSR J1550-5418 confirms that extremely high velocities need not result 
as a consequence of the formation of magnetars, and suggests that the velocity distribution 
of magnetars may not differ significantly from that of ordinary, rotation-powered neutron 
stars. 



The Long Baseline Array and the Parkes Observatory are part of the Australia Tele- 
scope, which is funded by the Commonwealth of Australia for operation as a National Facility 
managed by CSIRO. 
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Fig. 1. — Motion of PSR J1550-5418, with measured positions overlaid on the best fit, and 
shown as offsets from the best-fit position at MJD 54795.0 (Table [5]). (Left) Motion in 
Right Ascension vs time. (Right) Motion in declination vs time. 



